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Abstract. We demonstrate high-fidelity Zeeman qubit state detection in a single 
trapped *^S'r+ ion. Qubit readout is performed by shelving one of the qubit states to 
a metastable level using a narrow linewidth diode laser at 674 nm followed by state- 
selective fluorescence detection. The average fidelity reached for the readout of the 
qubit state is 0.9989(1). We then measure the fidelity of state tomography, averaged 
over all possible single-qubit states, which is 0.9979(2). We also fully characterize 
the detection process using quantum process tomography. This readout fidelity is 
compatible with recent estimates of the detection error-threshold required for fault- 
tolerant computation, whereas high-fidelity state tomography opens the way for high- 
precision quantum process tomography. 
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1. Introduction 

One of the basic requirements for implementing a physical qubit is the ability to 
faithfully measure its state |1]. Furthermore, qubit state detection must be performed 
with high fidelity to reach fault tolerant quantum computation. The exact detection 
error threshold required depends on many factors; however under quite general 
assumptions error values on the order of 10^^ — 10^"^ were estimated [2]. Another 
important use of high-fidelity state detection is precision process tomography for 
studying different quantum processes. 

Trapped ion qubits are a promising candidate system for physically realizing a 
quantum computer and serve as a convenient test-ground for studying fundamental 
quantum dynamics. Ion species used for this purpose typically have a single electron 
in their valence shell and the two qubit states are encoded in two energy levels of the 
valence electron. State detection methods used so far rely on state-selective fluorescence. 
Here photons are scattered from a laser that is resonant with a transition from one of the 
qubit states to a short-lived excited state, whereas transitions from the other qubit state 
are largely off-resonance. State inference is then based on the detected photon statistics 
[3]. To this end, qubit choices with a large energy separation are advantageous. 

For optical qubits, in which the two qubit states are separated by an optical 
transition, state detection fidelity as high as 0.9999 was demonstrated using state- 
selective florescence and accounting for photon arrival times [U |5]. Another detection 
scheme, using repetitive quantum nondemolition measurements on an ancila ion-qubit, 
was shown to give a fidelity of 0.9994 [6J. Optical qubits, however, have the disadvantage 
of an excited state lifetime on the order of one second, depending on the exact ion 
species used. Furthermore, the linewidth of even ultra-stable optical local oscillators, 
i.e., frequency stabilized lasers, would limit the dephasing time of an optical qubit. 

Qubits that are encoded into a pair of levels in the electronic ground state, in which 
the two qubit levels are split, either by the Zeeman effect or the hyperfine interaction, 
by radio-frequency transitions, have practically an infinite lifetime, as well as a very 
long coherence time [7]. Hyperfine qubit levels are typically separated by frequencies in 
the GHz range. State detection for this type of qubit can still be implemented by direct 
state-selective fluorescence, since the typical hnewidth of electric-dipole transitions is 
two orders of magnitude smaller. Here, off-resonance scattering normally limits state 
detection fidelity to below 0.995 [U [HJ [10]. The use of ancila qubits can, in principle, 
increase the detection fidelity of a hyperfine ion-qubit [11] . In ion species that have low- 
lying meta-stable levels, one of the qubit states can be shelved to a meta-stable level 
prior to detection. Here detection fidelity is similar to that of an optical qubit with an 
additional error introduced by the state shelving process. Using state selective optical 
pumping for shelving, a hyperfine ion-qubit measurement fidelity as high as 0.9977 was 
demonstrated. |1]. 

lon-qubits that are encoded into a pair of Zeeman split levels pose the hardest 
state-measurement challenge. This is because the frequency separation between the 
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qubit levels is typically in the MHz range and is comparable to the spectral linewidth 
of electric-dipole transitions used for state-selective fluorescence. Since direct state- 
selective fluorescence is impossible, Zeeman qubits can be readout only by state mapping 
onto a Hyperflne ancila qubit [6] or by shelving one of the qubit states onto a meta- 
stable level. However, the fldelity of state shelving in a Zeeman qubit using optical 
pumping is limited due to the strong coupling of the shelving light to the other qubit 
state [12]. Thus, the only way to state-selectively shelve a Zeeman qubit with high 
fldelity is by using a narrow-linewidth laser. Although there have been several reports 
on detecting a Zeeman qubit with shelving via a narrow-linewidth laser, with typical 
fldelities below 0.996 [IHl E], to our knowledge, there hasn't been a systematic study 
of the measurement error and limitations in this kind of qubit. 

Here we demonstrate the readout of a single-ion Zeeman qubit with a fldelity 
of 0.9989(1). Spin state detection is performed via electron shelving with a narrow- 
linewidth diode laser, followed by state-selective fluorescence. We analyze the different 
fundamental sources as well as technical sources of measurement error in detail. Further, 
we use the high-fldelity state detection of our qubit to demonstrate high-fldelity state 
tomography [15]. Here, we measure the fldelity of quantum state tomography averaged 
over all possible single qubit states and also fully characterize the detection process 
using quantum process tomography. 

2. Experimental Setup 

We trap a single ^^Sr~^ ion in a linear RF Paul trap. The trapping potential is well 
approximated as harmonic with a secular axial frequency of Uax = (271)1.09 MHz, and 
two nearly degenerate radial frequencies of Ur = (27r)2.5 MHz. A scheme of the relevant 
energy levels in ^^Sr~^ is shown in Fig. [H The | t) and | J,) qubit states are encoded 
in the 55*1/2,4.1/2 and 55*1/2-1/2 spin states respectively. A magnetic fleld of 0.477 mT 
splits the two qubit levels by a frequency, coq = (27r) 13.366 MHz. 

An oscillating magnetic fleld, perpendicular to the quantization axis, is generated 
by a current, oscillating at uo/{27r), through an electrode positioned 2 mm from the 
ion, resulting in coherent coupling between the two qubit levels. Zeeman qubit coherent 
rotations are therefore generated by oscillating current pulses. 

In the flrst part of state detection the electron is shelved from the | state to 
one of the Zeeman-split, levels using a narrow linewidth, 674 nm, external cavity 

diode laser. The diode laser is stabilized to a high flnesse (10^), ultra-low expansion 
glass, reference cavity. Residual frequency noise of the laser has three dominant spectral 
features. First, slow thermal drifts of the cavity result in drifts of the laser frequency 
with a typical magnitude of 10 Hz/S. The atomic resonance frequency is scanned every 
two minutes to correct for this drift. Second, intermediate frequency noise results in a 
laser line half-width of 70 Hz over several seconds, estimated in a Ramsey spectroscopy 
experiment [16]. The third spectral feature results from fast frequency noise that is 
unsuppressed, or even slightly enhanced, by our frequency servo system. This spectral 
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Figure 1: Energy levels scheme of the single valence electron of the ^^Sr'^ ion. The | '|") 
and I D qubit states are encoded in the two, Zeeman-split, spin 1/2 states of the 5*1/2 
ground level. Energy level lifetimes are written next to their spectroscopic notation. 
Laser light at 422 nm performs laser-cooling and state-selective fluorescence. Lasers at 
1092 nm and 1033 nm pump out population from the meta-stable -D3/2 and -D5/2 levels 
respectively. A 674 nm narrow linewidth diode-laser shelves the electron from the 
qubit levels to levels in the -D5/2 manifold. 



feature is often referred to as the servo bump. Our servo bumps are centered around 
700 kHz on both sides of the carrier, having a width of ~ 300 kHz , and contain an 
estimated < 5% of the optical power. A detailed discussion of our narrow linewidth 
diode laser system is given in [17]. In this experiment, the direction of 674 nm laser 
light propagation relative to the Zeeman splitting magnetic field direction, limits the 
allowed transitions between S1/2 and -D5/2 manifolds to Am = ±1 transitions. 

Electron shelving is followed by state-selective fluorescence detection on the 5*1/2 
Pi/2 transition at 422 nm. Photons that are scattered in a direction perpendicular to 
the 422 nm laser beam and the magnetic field are collected by an imaging system (N.A. 
= 0.31) and are detected with a Photo Multiplier Tube (PMT). State inference relies 
on the detected photon statistics. Qualitatively, a small number of detected photons 
implies that the ion qubit started in the | t) state and was shelved to the non-fluorescing 
1^5/2 level, whereas a large number of detected photons implies that the ion qubit started 
in the | l) state and therefore was not shelved and remained at the fluorescing 5*1/2 level. 

3. State Discrimination 

Following electron shelving, the number of photons n, detected by the PMT during a 
given detection time, tdet, is a random variable. This random variable is denoted by Uf, if 
the ion is in the fluorescing (bright) state 5*i/2, and if the ion is in the non-fluorescing 
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(dark) state Photon detection events, which occur when the ion is in the dark 

state are primarily due to scattering of the laser beam from trap surfaces. The fidelity 
of state discrimination is compromised by the overlap of the probability distribution 
functions (PDFs) of these two random variables. State inference can be performed by 
introducing a threshold value for the number of photons detected, nth- If the number of 
photons detected is greater (smaller) than this threshold, n > {<)nth, then we can infer 
that the ion is in the bright (dark) state. Given the probability distribution functions 
for nb and n^, the errors in detecting the bright and dark states are = Pb{n < nth) 
and erf = pd{n > nth) respectively. We want to find the parameters t^et and nth that 
minimize the mean error 

_ (^b + (^d _ Pb {n < nth) +Pd{n> nth) 

The detection fidelity is then given by F = 1 — e. If the lifetime of the -D5/2 level 
would have been infinite, the random variables nb and would follow two Poisson 
distributions. Given photon detection rates Rb and Rd in the bright and dark states 
respectively and a detection time, tdet, the means of these distributions would be given 
by ^b,d = Rb,dtdet- Here the longer the detection time, the smaller is the overlap between 
the two PDFs and therefore also the detection error. However, the finite lifetime, td^^^, 
of the -D5/2 level introduces a correction to the PDF for n^,, since there is a finite 
probability for the ion to decay during the detection. Upon decay, the photon detection 
rate becomes Rb- At a detection time much shorter than the lifetime, tdet <^ ''"D5/2; 
the PDFs for and are given by [18] 



Pb (n) = Poiss (n, nb) , (2a) 

Pd in) = f 1 - ^) Po^ss (n, n^) + ^ P (n„ n + 1) - F (n„ n + 1) ^ 
\ T J T nb-nd 

respectively. Here Poiss (n, n) denotes the probability of detecting n photons for a 
Poisson distribution with a mean n, and F (x, a) = Jq e~'^t°'~^dt is the incomplete 
gamma function. The two PDFs are denoted by the blue and red curves in Fig. [21 for 
our measured photon detection rates Rb = 73.5kHz and Rd = l.75kHz, a detection 
time of tdet = 285 fis and the known -D5/2 level lifetime, Tn^^^ = 390 ms [19] . 

As the detection time is increased, the overlap between the two functions initially 
decreases owing to the larger spacing between the two Poisson peaks, but it eventually 
increases owing to the growing tail of the dark distribution. Hence, an optimal detection 
time and a threshold number of photons exist such that the error in state discrimination 
is minimal. Figure [3a] shows a contour plot for e, as a function of the detection time tdet 
and the threshold on the number of photons nth, for the same Rb, Rd and td^^^ values 
used in Fig. [21 A minimal error of e ~ 2.9 ■ 10~^ is calculated at a detection time of 
tdet = 280//S and a nth = 5 threshold value for the number of detected photons. 

The actual distributions measured in the experiment are also affected by the state 
preparation error, and an error resulting from an imperfect shelving of the ion to the 
metastable state. State preparation errors, e^^inu and e^^inu for the | I) and | t) states, 
respectively, are given by the fraction of experiments in which the ion was initialized 
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number of photons 



Figure 2: Theoretical probability density functions for the number of detected photons 
in the bright (blue line) and dark (red line) states, given in Eqs. Eal and [26l using our 
measured photon detection rates Rf, = 73.5kHz and Rd = 1.75kHz, a detection time 
of tdet = 285 fis and the known level-lifetime, td^^^ = 390 ms. 



in the wrong state. The shelving error for the | t) state, e^^sheiving, is the probability 
that the ion remained at the 5*1/2 level after shelving was performed. For the | \.) state, 
^l,sheiving IS the probability that the ion was shelved to the level due to off- resonant 
light. Neglecting terms that are second order in the different errors, the resulting PDFs 
for the I I) and | t) states are given by [18], 

Pl (n) ~ (1 - ei^tot) Poiss (n, Ub) + ei^totPoiss (n, fid) (3a) 

Pt (n) ~ ^1 - e^^tot - Poiss (n, n^) + 

tdetT{nbn + l)-T{nd, n+1) r> ■ ^ - \ ^ot^ 
= = h e^^totPoiss [n, rib) ■ (3 b) 

T Ub-Ud 

Note that only the sum of the initialization error and the shelving error appears 
e^t-tot = (^in,init + ^in,sheiving- This prcvcuts US from distinguishing the state preparation 
error from the shelving error. The total mean detection error, 

- ^ e; + e^ ^ Pi {n < n^) +p^{n> n^) , . 

'2 2 ' ^ ' 

can be related to the error e resulting from the finite lifetime of the metastable level 
alone, obtained before 

-, ^l,tot + ^,tot 

e = e + ^ • (5) 

In particular, a minimum of e is obtained for the same values of detection time tdet and 
the threshold on the number of photons nth, as the minimum of e. 

4. Experimental Sequence 

In the experiment, two sets of data are taken. In each set the ion is prepared in one of 
the two Zeeman qubit states, | or | J,), and then state detection is performed. Each 
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Figure 3: State discrimination error as a function of the detection time tdet and the 
threshold on the number of photons Ufh- (a) Theoretical estimation obtained using 
Eqs. [2al [261 The measured photon detection rates of Rb = 73.5kHz Rd = 1.75kHz are 
assumed, as well as the known lifetime of the D5/2 level, t/j^^^ = 390 ms [19j. A 
minimal error of e = 2.9 ■ 10"^ is calculated at a detection time of tdet = 280/xs and a 
rith = 5 threshold value for the number of detected photons, (b) Experimental results. 
Here the mean initialization and shelving error of 8 ■ 10~^ was subtracted. A minimal 
error of e = 3(1) ■ 10~^ is measured at a detection time of tdet = 285/is and a threshold 
value for the number of photons of nth = 6. The theoretical model and our data are 
seen to be in relatively good agreement. 



set of data contains 3 ■ 10^ repetitions of the experiment, yielding statistical uncertainty 
of the estimated measurement error which is below 1 ■ 10"'^. 

In both sequences the ion is first Doppler cooled on the S1/2 — ?■ -P1/2 transition 
for 300 us, resulting in a mean axial harmonic oscillator number of n ^ 25. Second, 
sideband cooling is performed on the S'1/2,+1/2 — > D5/2,+3/2 narrow transition. To this 
end, the red sideband of the 5*1/2, +1/2 — > -D5/2,+3/2 transition is continuously excited for 
5 ms with the 674 nm laser, while the 1033 nm laser is left on to repump the population 
from the -D5/2 metastable level, and a polarized 422 nm light is left on to repump the 
population from the S'1/2,-1/2 state via the P1/2 manifold. Following sideband cooling 
the mean axial harmonic oscillator number is n = 0.3(2). 

Following cooling, the spin state is initialized. At this stage, it is highly probable 
that the qubit is already in the | state. To increase this probability further, we leave 
the polarized, 422 nm, optical pumping beam for an additional 50 fis. To initialize 
the qubit in the | 1) state, we use an 8 /is coherent qubit rotation to bring the electron 
from I 'I") to I 4,) . Optical pumping using the 674 nm laser follows, to further increase the 
initialization efficiency. Ten consecutive vr-pulses on the 5*1/2,^1/2 -D5/2,±i/2 transition, 
each followed by a 1033 nm repump pulse, pump the remaining population out of the 
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I 4,) or I t) state respectively. 

State detection begins with electron shelving. A 8.5 fis long, vr-pulse on the 
I 'I') — )■ -D5/2,+3/2 transition is applied. To increase shelving efficiency another, 14 /xs 
long, vr-pulse is applied on the | D^/2,-1/2 transition. Following shelving, an on- 

resonance 422 nm laser light is shined on the ion for 500 us, during which fluorescent 
photons are collected by the PMT and their time of arrival is recorded for further 
analysis. Then, any population shelved to the level is repumped back to the 

ground state using a, 100 /is long, 1033 nm laser pulse. At the end of the sequence, a 
red detuned 422 nm light Doppler-cools the ion until the next sequence begins. During 
the entire sequence the 1092 nm laser is left on to repump population from the -D3/2 
metastable level. 

5. State Detection Results 

Normalized histograms of the number of detected photons that were obtained in the 
two experiments for tdet = 285 fis, are plotted in Figs. Ha] and Hbl A maximum 
likelihood flt to the expected distribution functions given by Eqs. [3a| and [SJ] is denoted 
by the solid red line. The known detection time tdet = 285 fis and -D5/2 level lifetime 
''"^5/2 ~ 390 ms [19] are used whereas nd, ^i,tot and e^^tot are flt parameters. The sum 
of the initialization and shelving errors for the bright and dark states obtained from the 
flt are ei^tot = 6(1) ■ 10~^ and ei^^tot = 10(1) • 10~^ respectively. 

To determine the minimal detection error and the optimal parameters required 
to obtain it, we contour plot the measured mean error as a function of the detection 
time, tdet, and the threshold value, nth, in Fig- ES Here the mean error resulting from 
initialization and shelving = 8(1) ■ 10~^ is subtracted. As shown, the experimental 
plot reproduces the theoretical error plot shown in Fig. [3a] relatively well. In particular, 
the optimal parameters determined experimentally [tdet = 285/is, nth = 6) and the 
minimal error resulting from imperfect state discrimination e = 3(1) ■ 10~^ approach the 
estimated optimal parameters (tdet = 280/is, nth = 5), and detection error e = 2.9- 10"'*. 

6. Error Budget 

The error shown in Figures l3a] and [3b] is fundamental and does not result from technical 
imperfections. As shown by Myerson et. al. [1], this error can be somewhat reduced by 
accounting for the photon times of arrival 0. Shelving and initialization errors however, 
result from technical imperfections. In the following discussion we try to point out 
the dominant sources for this errors by calculating an order of magnitude estimates for 
the contribution of different possible error sources. A summary of this error budget is 
presented in Table [TJ 

I By similarly analyzing photon arrival times we were indeed able to lower this error to 2(1) • 10~^ 
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Figure 4: Experimentally obtained PDFs for a detection time of tdet = 285 fis. The ion 
is initialized in the | \.) state in Fig. (a) and in the | state in Fig. (b). A maximum 
likelihood fit to the expected distribution functions, given by Eqs. [3a| and |3fol are 
shown (solid red line). Total initialization and shelving errors for the | \.) and | t) 
states obtained from the fit are = 6(1) • 10~^ and = 10(1) ■ 10~^ respectively. 



6.1. Initialization Error 

The first stage of initialization consists of optical pumping to the | t) state via the P1/2 
manifold. Ideally, the optical pumping beam has to match a perfect polarization. 
Otherwise, the matrix elements for the Si/2,+1/2 -Pi/2,±i/2 transitions will not null, 
and some population will be pumped out of the | t) state. To this end, the wave-vector, 

— * 

k, of the beam has to be exactly parallel to the external magnetic field and a perfect 
circular polarization is required. We found the fidelity of this initialization step in our 
setup to be limited to less than 0.999 due to stress-induced birefringence in the fused 
silica vacuum chamber optical ports. 

To increase initialization efficiency, we use the narrow linewidth 674 nm laser 
to optically pump the remaining population from the | l) state to the | t) state via 
the -D5/2,+i/2 level as described above. This initialization procedure is limited due to 
the finite off-resonance excitation probability on the S'1/2,+1/2 — > D5/2,+3/2 transition, 
resulting in a finite probability for the electron to be in the | \.) state at the end of 
the process. Note that there are two different contributions to off-resonant excitation. 
The first is coherent with respect to the resonant light component and is due to the 
pulse finite time. The second is incoherent and is due to the servo bumps. In steady 
state, i.e., after a sufficient number of pulses, this probability is given by the ratio of the 
I t) — 7- I 4,) and I i) — I t) transfer probabilities. These are, in turn, estimated based 
on the off-resonance excitation rate, the polarization of the 1033 nm repump laser, and 
the different decay probabilities from the P3/2 excited states to either the | t) or I i) 
states. The incoherent off-resonance excitation rate is directly measured and is seen 
to be very sensitive to the specific laser and servo-loop working parameters, such as 
current or gain; we have observed this rate change by a factor of 2 when working with 
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slightly different parameters. The coherent off-resonance excitation rate is estimated 
based on the measured Rabi frequencies and detuning. We thus estimate the | t") 
state initiahzation error to be ~ 1 ■ lO^'' due to coherent off-resonance excitation and 
~ 0.5 ■ 10^'^ due to incoherent off-resonance excitation. 

When initiahzing to the | I) state, we initially perform optical pumping to the 
I t) state and apply an RF TT-pulse, which transfers the electron to the | 4-) state with 
a fidelity higher than 0.999. This procedure is followed by 10 similar pulses on the 
'S'i/2,+1/2 ~^ D5/2-1/2 transition, which result in | 4,) state initialization with the same 
error as for the | t) state. 

Note that another possible error in the initialization process could arise from leakage 
of 422nm laser light, resulting in mixing of the two qubit states. This was, however, 
found not to be a problem in our system by measuring the lifetime of each of the qubit 
states. 

6.2. Electron Shelving Error 

The error in the shelving process is different for the two qubit states and is therefore 
analyzed separately. 

I D state shelving error. For the | I) state, the shelving error is the probability that 
the ion was shelved to the -D5/2 level manifold due to off-resonant light. The nearest 
transition from | 1) to the -D5/2 level is 5.35 MHz detuned from the laser carrier, and 
the coupling to the motional sidebands is suppressed by the Lamb-Dicke parameter 
{t] = 0.05 for the longitudinal motion). Here, the contribution of coherent off-resonance 
excitations during the two shelving pulses to the error is estimated to be ~ 2 • 10"^, 
while the error due to incoherent excitations is ~ 1 • 10~^. 

I t) state shelving error. For the | t) state, the shelving error is the probability that 
the ion remains in the S1/2 manifold after shelving was performed. A number of factors 
contribute to this error. 

The shelving transition Rabi frequency depends on the ions' motion through the 
Debye- Waller factor. The temperature of the ion therefore has a large effect on the 
shelving error. The ion is initially Doppler cooled to a mean axial harmonic oscillator 
number of n ^ 25. This leads to an error of ~ 0.01 in a single shelving 7r-pulse. It is 
important to note that in the second shelving pulse the error is increased to ~ 0.1 due 
to the change in the level occupation distribution induced by the first pulse; following 
the first pulse, motional states with a smaller Debye- Waller factor are preferentially 
left at the S'1/2 level. To reduce this error, Doppler coohng of the ion is followed by 
sideband cooling on the <S'i/2,+i/2 — > -D5/2,+3/2 transition. After five ms cooling, a mean 
axial harmonic oscillator number of n = 0.3(2), inferred from the imbalance between 
the motional sidebands, is reached. Yet, the mean harmonic oscillator level does not 
provide the full motional population distribution, which is required to estimate the 
shelving error. A thermal distribution with n = 0.3 would lead to a shelving error 
equal to ~ 1 • 10~^ in the first pulse and ~ 1 ■ 10^"^ in the second pulse. However, 
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the distribution of high harmonic oscillator levels, following sideband cooling, is poorly 
described by a thermal distribution. This is because motional states outside the Lamb- 
Dicke regime have a small Debye- Waller factor on the cooling transition and therefore 
are not efficiently cooled. The exact dynamics of the sideband cooling process and 
the resulting energy level distribution are not calculated here. In our estimates we use 
the error values assuming a thermal distribution. This is probably an overly optimistic 
estimate and could be the source of the discrepancy between our evaluated and measured 
errors for the | state. 

Motion along the radial direction of the trap has a much weaker effect on shelving 
infidelity, because of the stronger confinement along this direction. However, while 
sideband cooling is performed for the axial motion, the radial motion is only Doppler 
cooled and thus can not be neglected. We estimate the infidelity due to radial motion 
to be ~ 3 ■ 10""^ in the first pulse and ~ 1 ■ 10~^ in the second pulse. 

Intensity fiuctuations of the 674 nm laser are reduced to a negligible level using an 
intensity noise eater. An effective residual intensity noise resulting from beam pointing 
fiuctuations on the ion is estimated to give an error of ~ 3 ■ 10~^ in a single shelving 
pulse. 

Magnetic field fiuctuations are studied in a Ramsey experiment on the | t) I i) 
transition, where a coherence time of 500/is is measured. The resulting errors are 
estimated by numerical solution of the optical Bloch equations with a corresponding 
coherence decay term. Taking into account the different magnetic field susceptibilities 
for the different transitions, as well as the different durations for the two shelving pulses, 
we estimate the errors to be ~ 1.5-10"^ and ~ 4.5-10"^ for the two vr-pulses respectively. 

A different contribution to the shelving error results from laser frequency noise. 
As detailed above, this noise has few spectral contributions. Noises that are slow 
compared with the experiment time result in a typical laser frequency drift of two 
kHz between consecutive scans of the atomic line, giving a shelving error of ~ 3 ■ 10~^ 
per pulse. The contribution of the laser linewidth to the decoherence rate is found in a 
Ramsey experiment on the 5*1/2, +1/2 — ^ -D5/2,+3/2, yielding a coherence time of 700 /is. 
After quadratically subtracting the contribution of magnetic field fiuctuations to this 
decoherence rate, we estimate the error owing to the linewidth of the laser alone to be 
~ 2 • 10~^ in a single shelving pulse. 

The finite decay probability from the -D5/2 level manifold during the shelving process 
also contributes to the | t) state shelving error. This probability has two contributions. 
The first is due to the finite lifetime of the metastable level rn^,^ = 390 ms [T9l. After 
the first shelving pulse, the electron is found at the level with high probability 
and therefore it decays with rate — - — during the second shelving pulse. During the 

first shelving pulse this error is half as large, because, on average, only half of the 
population is in the -D5/2 manifold. Using pulse durations of 8.5 fis and 14 /is we 
estimate this effect to give a contribution of ~ 0.5 ■ 10~^ to the shelving error. The 
second contribution is due to off-resonant coupling of light to motional sidebands. This 
coupling results in population transfer from the manifold back to the | t) state. 
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The largest such contribution is due to the incoherent light in the servo bump which has 
a significant spectral overlap with the axial-motion sidebands, located 1.1 MHz away 
from the carrier. The estimated contribution to the error here is ~ 0.5 ■ 10""^ for the two 
shelving pulses combined. 

In addition, leakage of lOSSnm laser light during detection can shorten the lifetime 
of the meta-stable state and thus reduce detection fidelity. In this experiment we 
find the lifetime of the -D5/2 level to be consistent with [19], indicating that leakage of 
repump light has no significant effect. 

Summing up all the different contributions for the | '|") state shelving error after 
both TT-pulses we get an error estimate equal to ~ 2.5 ■ 10~^. 
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Shelving Error Sources 
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2 


Incoherent off-resonance excitation 




1 


Off-resonance sideband excitation 


0.5 




-D5/2 level decay during shelving 


0.5 




Axial motion in the trap 


0.1 


1 




Radial motion in the trap 


3 


10 




Beam pointing noise 


30 


30 




Frequency drift 


30 


30 




Fast laser linewidth 


20 


20 




Magnetic field noise 


15 


45 




Shelving Total Error 


2.5 


3 


Overall estimated error 


4 


4.5 


Overall measured error 


10(1) 


6(1) 



Table 1: Shelving and initialization error budget for the two different qubit states. 
Errors are given in units of 10~^. The initialization error results from off-resonant 
coupling to the other qubit state. The different error sources in the shelving operation 
are detailed. The total shelving error for the | t) state is obtained by summing all 
shelving error sources in each of the pulses, excluding the contributions of different 
off- resonance excitations and -D5/2 level decay. After the two sums are multiplied, the 
excluded errors are added. 
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6.3. Error Summary 

A summary of the estimated initialization and shelving errors and their sources is 
presented in Table [1] in units of 10"'^. The total estimated errors for the | ^ind 
I D states are ~ 4 ■ 10~^ and ~ 4.5 ■ 10~^ respectively, whereas the measured errors of 
10(1) ■ 10"'* and 6(1) ■ 10~^ respectively are roughly twice as large. 

As the error budget suggests, one significant source for initialization and shelving 
errors is off-resonance excitation. This is not surprising considering that here one needs 
to spectrally distinguish between two states that differ by a 13 MHz on top of a 445 THz 
transition. There are two main contributions to off-resonance excitation. The first is 
coherent and is due to the pulse finite time, while the second is incoherent and is driven 
by the laser servo bumps. There are several possible avenues for reducing these errors. 
Increasing the magnetic field and thus the separation between the qubit states will 
reduce the off-resonance excitation probability. Decreasing the shelving laser intensity 
will reduce both coherent and incoherent off-resonance excitation errors. Slowing down 
the shelving transition Rabi frequency, thereby increasing the pulse length, will reduce 
the coherent off-resonance excitation error in proportion to the laser intensity. Second, 
since the pulse time is proportional to the square root of the laser intensity whereas 
the incoherent off-resonance excitation rate is proportional to the laser intensity, the 
incoherent off-resonance excitation error will reduce as the intensity square root. Note 
that this error reduction will come at the expense of larger error contributions due to 
slower laser frequency noise (drift and linewidth) and magnetic field noise. The shelving 
laser servo bumps can be reduced by engineering a faster servo system, by using a laser 
that has a narrower noise bandwidth to begin with (diode laser systems are notorious 
for their broad frequency noise spectrum), or by spectral filtering. 

Here, incoherent off-resonance excitation errors are below 10~^. However, it is 
worth noting that this error source becomes much larger when carrying out operations 
that are off-resonance from the carrier, where the servo bump and the transition carrier 
have a significant overlap. We observe a large error caused by the servo bump when 
performing Rapid Adiabatic Passage (RAP) on the 5*1/2 -D5/2 transition [13], as well 
as when driving the motional sideband. 

As mentioned above, a possible explanation for the discrepancy between the 
estimated and measured | state errors is a deviation of the ions' harmonic oscillator 
energy level distribution from a thermal distribution. Such a deviation is not surprising 
since in instances in which the ion motion is high, both the sideband cooling mechanism 
and the shelving transition are likely to fail. As an example, starting from a thermal 
distribution after Doppler cooling, failure to sideband cool all the instances where the 
ion motion is in n > 90 will produce a larger shelving error than that we measure. 
Reducing this error contribution will require better cooling of the tail of the ion energy 
distribution, via e.g. second-sideband cooling. 
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7. High-Fidelity Quantum State Tomography 

lon-qubit state measurement is an important tool for executing many quantum 
algorithms [20l|2Tl[22] and for studying different quantum processes via state and process 
tomography [23l [24l |25]. In all these cases the states that are being measured do not 
necessarily coincide with the electronic eigen-basis (the two Zeeman states in the case 
of a Zeeman qubit). However, all these states are related to the electronic eigen-basis 
via single qubit rotations. Typically, measurement fidelity has been characterized as 
the average fidelity of state detection for the two electronic eigen-states (which are 
also the measurement eigen-basis). The fidelity of single-qubit rotations was separately 
benchmarked [26]. However, the error introduced by rotations is not uniformly spread 
over all possible input states. To this end, the experimental estimate of state tomography 
fidelity, averaged over all possible single ion-qubit states, is beneficial. 

Qubit state tomography is represented by a quantum map e. The fidelity of state 
tomography of a given pure input state pj = |\E'j)(\E'j| is therefore the fidelity of this 
state and the reconstructed output state e (pj), 

F = Tr{e{p,)p,). (6) 

The output state e (pj) is reconstructed by 

Here p^i Py and are the measured projections of pj on the x, y and z axes 
correspondingly and crj are the Pauli spin operators [I5]. Note that this definition 
of the fidelity is in agreement with Eq. [T]for the special cases of |\E'j) = | |), | f). 

Here we are interested in determining the fidelity of state tomography averaged over 
all possible input states. A value for the fidelity, averaged over all possible qubit states, 
can be obtained by calculating an algebraic average of the fidelities of the six pure input 
states I- ^) = \l),\ + z) = \ t), lix) = (| ;) ± | t))/\/2 and | ± ?/) = (| ;)±z| t))/\/2 
[28] . These measurements require the ability to initialize the qubit in different states, as 
well as perform measurements in different bases. Both are achieved via qubit rotations. 

The experimental sequence performed is similar to the one previously discussed. 
Briefly, the ion is flrst ground-state cooled. Then state preparation is performed; 
the qubit is flrst prepared in either the \ + z) or \ — z) state and when required 
is rotated to initialize the | ± x),\ ± ?/) states. Projection measurement consists of 
mapping the measurement basis onto the | ± z) basis, once again using qubit rotations, 
followed by shelving and state-selective fluorescence detection. In principle a total of 
18 measurements, three for each input state, are needed. However since there is no 
phase information in the detection process, only nine measurements are not redundant 
and were therefore performed. The measurement uncertainties are determined by the 
amount of collected statistics. For the projections on an axis that is parallel to the input 
state, 2 ■ 10^ repetitions were performed, yielding an error of 2 ■ 10~^. For projections 
on an axis that is orthogonal to the input state 1 ■ 10^ repetitions were performed which 
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Table 2: Summary of the state tomography results. The uncertainties in the values of 
the projections on axis parallel to the initial state are 2 ■ 10~^. For projections on the 
orthogonal axis the uncertainties are 1 • 10~^. The uncertainties for the calculated 
fidelities are presented in the table. 



yield an uncertainty of 1 ■ 10~^. A summary of the different measured projections and 
the calculated fidelities for different input states are presented in Table [2l 

The I ± z) measurement fidelities here are somewhat lower than the best effort 
presented above. This might be due to slightly non-optimal detection parameters 
The measurement fidelities for the | ± x), | ± y) states are generally lower than the 
I ± z) measurement fidelities due to larger state initialization and detection errors 
caused by imperfect qubit rotations. The non-zero projections in measurement basis 
orthogonal to the initialization axis are due to systematic errors in the preparation and 
the measurement sequence (e.g. slightly incorrect pulse durations). Note, however, 
that all the measured state projections are within two standard deviations from the 
expected value [Jj]. Based on the measurements performed, we calculate the averaged 
state tomography fidelity over the entire Bloch sphere to be F = 0.9979(2). 

A more complete characterization of the detection process is achieved by performing 
full process tomography [15]. We use the chi matrix, x, representation to characterize 
the completely positive map representing the detection process, 

4 

e (Pj) = J2 XmnEmpjEl (8) 
m,n=l 

Here, the fixed set {i?m}m=i = {I,o'x,'io'y,o'z} forms a basis for single qubit quantum 
maps. Ideally, the detection process would be represented by the chi matrix 
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(9) 



§ The optimal tdet and nth values depend on the photon detection rates, Rb and Rd, which can drift 
due to small changes in the laser parameters. 

II In fact, taking into account the imbalance in the fidelities of ±Z , the projection on an orthogonal 
axis is expected to be —2 • 10"'' 



High-Fidelity State Detection and Tomography of a Single Ion Zeeman Qubit 



16 



corresponding to the Identity operation. 

To determine the chi matrix of a single qubit map experimentally, it is enough to 
measure the output density matrices for the following set of four, linearly independent, 
input density matrices, {| + z){-'rz\, \ — z){—z\, \ + x){-\-x\, \ + y){+y\}. An exphcit 
formula can then be obtained for the chi matrix values [I5]. The non-zero orthogonal 
state projections, originating from systematic and statistical errors, result in a non- 
physical (non-positive) reconstructed chi matrix. To obtain a meaningful physical chi 
matrix, we therefore null all orthogonal state projections 0. The absolute values of 
the reconstructed chi matrix entries are plotted in Fig. [5a] on a logarithmic scale. As 
expected, xii? which represents the identity operation, is three orders of magnitude 
larger than any other entry. Other diagonal entries are an order of magnitude larger 
than off-diagonal entries, implying that qubit depolarization is the main error channel. 
The increased measurement error of superposition states, due to rotation errors, is 
manifested in the slightly larger X44 (dephasing channel) as compared with X22 and X33 
(spin- flip channels). Off-diagonal elements are due to the small imbalance between the 
measurement fidelity of the | + z) and the | — z) states. 

Using the obtained chi matrix, one can calculate the output of any input state, 
and hence also the fidelity. Figure [5b] depicts the detection error over the entire 
Bloch sphere. The polar angle, 6, and the azimuthal angle, 0, define the pure state 
|\E') = cos (^1^ I t) + 6**^5^(^1^ I l), whereas the radial distance from the origin in 
this direction indicates the error for this state. As expected from qubit depolarization 
an almost spherical spheroid is obtained. The resulting spheroid is also somewhat 
oblate due to the added error in qubit rotations and is minutely dilated in the positive 
hemisphere direction, owing to the slight imbalance between the measurement fidelity 
of the I + z) and | — z) states. Using the reconstructed chi matrix, we find the process 
fidelity to be Fp^oc = Tr{xideaiXproc) = 0.997(1). 

8. Summary 

In conclusion, we demonstrate high-fidelity state detection of a qubit encoded in the 
Zeeman splitting of the ground state of a single ^^Sr~^ ion. The limitations of the 
best effort readout fidelity of 0.9989(1) are analyzed in detail. Our estimates of the 
contribution of imperfect state preparation and shelving to the measured error is 
~ 8 ■ 10~^. The remaining part, ~ 3 ■ 10~^ , results from limited state discrimination 
efficiency owing to the finite lifetime of the metastable level. This fraction of the error is 
intrinsic to our state detection scheme, and can be somewhat decreased if the information 
on photon detection times is used [1]. 

^ A different approach would have been to perform maximum hkehhood estimate of a physical operation 
[251 [27]. However since we believe we have identified the origin of the operation non-positiveness with the 
small finite projection in the orthogonal direction, the better approach, in our opinion, is to null those. 
Two more notes: with better statistics the magnitude of orthogonal projections can be increasingly 
lowered by adjusting the rotation pulse time. Second, by artificially nulling these projections the 
estimate of the error becomes larger and is therefore pessimistic. 
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Figure 5: Process tomography results for the detection process, (a) Absolute values of 
the reconstructed chi matrix entries, (b) Detection error for every pure input state 
reproduced using the reconstructed chi matrix. Here, the error in detecting the pure 
state, associated with given azimuthal and polar angles on the Bloch sphere, is 
represented by the radial distance from the origin of the shown surface in that 
direction. The resulting spheroid is slightly oblate due to an increased error in the 
equatorial direction. This additional error is due to imperfections in the RF pulses 
used to initialize qubit superpositions. 

Since one of the applications of high-fidelity state detection is performing high 
accuracy quantum process tomography we measured the averaged state tomography 
fidelity over the entire Bloch sphere, which is 0.9979(2). We also performed quantum 
process tomography for the detection process, and found the process fidelity to be 
Fproc = 0.997(1). This fidelity can be further increased if higher fidehty qubit rotations 
are used. 

Our measured state detection fidelity is compatible with recent estimates of 
the required fault-tolerance error threshold and can be used in the future to 
implement effective quantum error correction protocols. In addition, highly accurate 
quantum process tomography, important for studying basic quantum processes, can be 
implemented. 

We gratefully acknowledge supported by the ISF Morasha program, the Crown 
Photonics Center and the Minerva Foundation. 
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